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Abstract. In this paper we study the existence of a unique solution to a general class of 
Young delay differential equations driven by a Holder continuous function with parameter 
greater that 1/2 via the Young integration setting. Then some estimates of the solution 
are obtained, which allow to show that the solution of a delay differential equation driven 
by a fractional Brownian motion (ffim) with Hurst parameter H > 1/2 has a C°°-density. 
To this purpose, we use Malliavin calculus based on the Frechet differentiability in the 
directions of the reproducing kernel Hilbert space associated with ffim. 



1. Introduction 

The recent progresses in the analysis of differential equations driven by a fractional 
Brownian motion, using either the complete formalism of the rough path analysis [3j [TOl 
[18], or the simpler Young integration setting [25], [33], allow to study some of the basic 
properties of the processes defined as solutions to rough or fractional equations. This 
global program has already been started as far as moments estimates [13], large deviations 
[T6] . or properties of the law [21 [21] are concerned. It is also natural to consider some 
of the natural generalizations of diffusion processes, arising in physical applications, and 
see if these equations have a counterpart in the fractional Brownian setting. Some partial 
developments in this direction concern pathwise type PDEs, such as heat [71 [TT1 [T2l 130] . 
wave [28] or Navier-Stokes [4] equations, as well as Volterra type systems [6]. As we 
shall see, the current paper is part of this second kind of project, and we shall deal with 
stochastic delay equations driven by a fractional Brownian motion with Hurst parameter 
H > 1/2. 

Indeed, we shall consider in this article an equation of the form: 

dy t = f(Z?)dB t + b{Z v t )dt, t G [0,T], (1) 

where B is a <i-dimensional fractional Brownian motion with Hurst parameter H > 1/2, 
/ : Cj([-h,0];R n ) -> R nxd and b : Cj([-h,0];R n ) -> R n satisfy some suitable regularity 
conditions, Cj designates the space of 7-Holder continuous functions of one variable (see 
Section |2~TT1 below) and Z\ : [— h, 0] — > IR n is defined by Zf(s) = yt+ s - In the previ- 
ous equation, we also assume that an initial condition £ G C\ is given on the interval 
[— h,0]. Notice that equation (CQ) is a slight extension of the typical delay equation which 
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is obtained for some functions / and b of the following form: 




o 



Vt+ev 



(d9) , (2) 



for a regular enough function a, and a finite measure v on [— h,0]. This special case of 
interest will be treated in detail in the sequel. Our considerations also include a function 
/ defined by f(2f) = cr(Z t y (— ui), . . . , 2%(— for a given k > 1, < u\ < . . . < Uk < h 
and a smooth enough function a : M. nxk —>■ M. nxd . 

The kind of delay stochastic differential system described by ([]]) is widely studied when 
driven by a standard Brownian motion (see [20] for a nice survey), but the results in 
the fractional Brownian case are scarce: we are only aware of [8] for the case H > 1/2 
and f(Z y ) = a(Z y (—r)), < r < h, and the further investigation [9] which establishes 
a continuity result in terms of the delay r. As far as the rough case is concerned, an 
existence and uniqueness result is given in [22j for a Hurst parameter H > 1/3, and [3T] 
extends this result to H > 1/4. The current article can be thus seen as a step in the study 
of processes defined as the solution to fractional delay differential systems, and we shall 
investigate the behavior of the density of the M n -valued random variable y t for a fixed 
t £ (0, T], where y is the solution to (DQ). More specifically, we shall prove the following 
theorem, which can be seen as the main result of the article: 

Theorem 1.1. Consider an equation of the form (T7]j for an initial condition £ lying in 
the space CJ([— h, 0]; W 1 ). Assume b = 0, and that f is of the form (HJ) for a given finite 
measure v on [— /i,0] and a : 1R" — > M. nxd a four times differentiate bounded function with 
bounded derivatives, satisfying the non-degeneracy condition 



Suppose moreover that H > H , where H = (7 + vTf)/16 ~ 0.6951. Let t £ (0,T] 
be an arbitrary time, and y be the unique solution to (T7]j in ([0, T]; M n ), for a given 
1/2 < k < H. Then the law of y t is absolutely continuous with respect to Lebesgue 
measure in M. n , and its density is a C°° -function. 

Notice that this kind of result, which has its own interest as a natural step in the study 
of processes defined by delay systems, is also a useful result when one wants to evaluate 
the convergence of approximation schemes in the fractional Brownian context. We plan 
to report on this possibility in a subsequent communication. The reader may also wonder 
about our restriction H > H above. It will become clear from Remark 13.151 that this 
assumption is due to the fact that we consider a delay which depends continuously on 
the past. For a discrete type delay of the form a(y t ,y t ^ n , . . . ,y t - r ), with q > 1 and 
r± < ■ • • < r q < h, we shall see at Remark 14.71 that one can show the smoothness of the 
density up to H > 1/2, as for ordinary differential equations. Finally, the case 6 = has 
been considered here for sake of simplicity, but the extension of our result to a non trivial 
drift is just a matter of easy additional computations. 

Let us say a few words about the strategy we shall follow in order to get our Theorem ll.il 
First of all, as mentioned before, there are not too many results about delay systems 
governed by a fractional Brownian motion. In particular, equation (pQ) has never been 
considered (to the best of our knowledge) with such a general delay dependence. We 
shall thus first show how to define and solve this differential system, by means of a slight 



for all ?7i, ?7 2 G M". 
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variation of the Young integration theory (called algebraic integration), introduced in [10] 
and also explained in [21]. This setting allows to solve equations like (CQ) in Holder spaces 
thanks to contraction arguments, in a rather classical way, which will be explained at 
Section 13.11 In fact, observe that our resolution will be entirely pathwise, and we shall 
deal with a general equation of the form 

dy t = f(Z?)dx t + b(Z?)dt, te[0,T], (3) 

for a given path x G ([0, T}; M. d ) with 7 > 1/2, where the integral with respect to x has 
to be understood in the Young sense |32j. Furthermore, in equations like ([3]), the drift 
term b(Z y ) is usually harmless, but induces some cumbersome notations. Thus, for sake 
of simplicity, we shall rather deal in the sequel with a reduced delay equation of the type: 

yt = a+ f f(Z y s )dx s , te[0,T}. 
Jo 

Once this last equation is properly defined and solved, the differentiability of the solution 
y t in the Malliavin calculus sense will be obtained in a pathwise manner, similarly to the 
case treated in [26]. Finally, the smoothness Theorem 11.11 will be obtained mainly by 
bounding the moments of the Malliavin derivatives of y. This will be achieved thanks to 
a careful analysis and some a priori estimates for equation (JT|). 

Here is how our article is structured: Section [2] is devoted to recall some basic facts 
about Young integration. We solve, estimate and differentiate a general class of delay 
equations driven by a Holder noise at Section [3j Then at Section 0] we apply those general 
results to fBm and prove our main Theorem 11.11 

2. Algebraic Young integration 

The Young integration can be introduced in several ways (convergence of Riemann sums, 
fractional calculus setting [33]). We have chosen here to follow the algebraic approach 
introduced in [10] and developed e.g. in [12 [21], since this formalism will help us later in 
our analysis. 

2.1. Increments. Let us begin with the basic algebraic structures which will allow us to 
define a pathwise integral with respect to irregular functions: first of all, for an arbitrary 
real number T > 0, a topological vector space V and an integer k > 1 we denote by C k (V) 
(or by Cfc([0, T]; V)) the set of continuous functions g : [0,T] fc — > V such that gt Y ---t k = 
whenever U = t i+ \ for some i < k — 1. Such a function will be called a (k — 1) -increment, 
and we will set C*(V) = Ufc>iCfc(^). An important elementary operator is 5, which is 
defined as follows on Ck(V): 

k+l 

5 : C k {V) - C k+1 (V), (Sg) tl ... tk+1 = If '//,,.,,.,,.,• (4) 

i=l 

where U means that this particular argument is omitted. A fundamental property of 5, 
which is easily verified, is that 55 = 0, where 55 is considered as an operator from CkiV) 
to C k+2 {V). We will denote ZC k (V) = C k (V) n Ker<5 and BC k {V) = C k (V) n ImS. 

Some simple examples of actions of 5, which will be the ones we will really use through- 
out the paper, are obtained by letting g G Ci(V) and h G C2CV). Then, for any 
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s,u,t G [0, T], we have 

(5g) st = g t - g s , and (5h) sut = h st - h su - h ut . (5) 

Furthermore, it is easily checked that ZCk(V) = BCk(V) for any k > 1. In particular, the 
following basic property holds: 

Lemma 2.1. Let k > 1 and ft G i?Cfc + i(V). Then there exists a (non unique) f G Ck(V) 
such that h — Sf. 

Observe that Lemma [2J] implies that all the elements h G C 2 (V) such that 5h = 
can be written as h = Sf for some (non unique) / G Ci(V). Thus we get a heuristic 
interpretation of S\c 2 (y)'- it measures how much a given 1-increment is far from being an 
exact increment of a function, i.e., a finite difference. 

Remark 2.2. Here is a first elementary but important link between these algebraic struc- 
tures and integration theory: let / and g be two smooth real valued function on [0, T]. 
Define then I G C 2 (V) by 




for s,te[0,T}. 



Then, some trivial computations show that 

(SL) sut = [g u - g s ][f t - f u ] = (Sf)ut(Sg)su- 

This is a helpful property of the operator S: it transforms iterated integrals into products 
of increments, and we will be able to take advantage of both regularities of / and g in 
these products of the form Sf 5g. 

For sake of simplicity, let us specialize now our setting to the case V = M. m for an 
arbitrary m > 1. Notice that our future discussions will mainly rely on ^-increments with 
k < 2, for which we will use some analytical assumptions. Namely, we measure the size 
of these increments by Holder norms defined in the following way: for < a\ < a 2 < T 
and / G C 2 ([ai, a 2 ]; V), let 

H/IL^aa] = SU P |, ^ ■„ > and C 2 M ([ai, a 2 ]; V) — •[/ G C 2 (V); Wf]]^,^] < °°} ■ 

r,t€[ai,a 2 ] I 1 r \ 

Obviously, the usual Holder spaces Cf ([ai, a 2 ]; V) will be determined in the following way: 
for a continuous function g G Ci([ai, a 2 ]; V), we simply set 

IML[ai,aa] = ||^IU[ai,oa]> ( 6 ) 

and we will say that g G Cf ([ai, 02]", V) iff \\g || ^,[01,02] i s finite. Notice that || ■ H^jai.oa] is 
only a semi-norm on Cf ([ai, a 2 ]; V), but we will generally work on spaces of the type 

c v, au a 2 ( v ) = {9 ■ [ai,a 2 ] -> V; # 01 = «, Hfl-IU,^,^] < 00} , (7) 
for a given v G V, or 

:= {C e Cf(K - h,a 2 ];R d y,C = Q on [a, - h, ai ]}, (8) 

where < a\ < a 2 and g G Cf ([ai — ft, ai]; These last two spaces are complete metric 
spaces with the distance d^(g,f) = \\g — /|| M . More specifically, the metric we shall use 
on the space Cg aia2 (M. d ) is: 

<^/x,ai,02 {.9i f) = Wd ~ f\\n,[ai— h,a,2\- 
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In some cases we will only write C%(V) instead of C£([ax, a 2 ]; V) when this does not lead 
to an ambiguity in the domain of definition of the functions under consideration. For 
h G Cs([ai, a 2 }; V) set in the same way 

WHi^M = S ^P \ u JsMt-u\P (9) 

s,u,te[ai,a 2 ] \ u b \ r u \ 



inf I J] ||/iilU,/i-«; ^ = ^' < Pi < ^ J 



where the last infimum is taken over all sequences {hi G C 3 (V)} such that h = 

and for all choices of the numbers G (0,/i). Then || • || M is easily seen to be a norm on 

C 3 ([ai, a 2 ]; V), and we set 

C%{[a u a 2 ]-V) := {h G C 3 ([a 1 ,a 2 ]]V); ||% < 00}. 

Eventually, let C 3 + ([ax,a 2 ];V) = U^iCg ([ai, 02]", V), and remark that the same kind of 
norms can be considered on the spaces ZC 3 ([ai, a 2 ]; V), leading to the definition of some 
spaces ZC% ([ai, a 2 ]; V) and ZC\ + ([a\, a 2 ]; V). 

With these notations in mind, the crucial point in our approach to pathwise integration 
of irregular processes is that, under mild smoothness conditions, the operator 5 can be 
inverted. This inverse is called A, and is defined in the following proposition, whose proof 
can be found in [TO] . 

Proposition 2.3. Let < a± < a 2 < T. Then there exists a unique linear map A : 
ZCl + ([ai, a 2 \; V) — > C\ + ([ai, a 2 ]; V) such that 

SA = Id zc i+ {[aita2] . vy 

In other words, for any h G Cg + ( [01,02]; V) such that 5h = there exists a unique g = 
A(h) G C2 + ([ai, 02]; V) such that 5g = h. Furthermore, for any [i > 1, the map A is 
continuous from ZC 3 ([ai, a 2 }; V) to C 2 ([a\, a 2 }; V) and we have 

\\Ah\\ heZC^[ ai ,a 2 ];V). (10) 

Moreover, the operator A can be related to the limit of some Riemann sums, which 
gives a second link (after Remark 12.21) between the previous algebraic developments and 
some kind of generalized integration. 

Corollary 2.4. For any 1-increment g G C 2 (V) such that 5g G C\ + , set 5f = (Id— A5)g. 
Then 

n-l 

(5f) st = lim 

\i-i-st — HJ . 

1=0 

where the limit is over any partition U st = {t = s, . . . ,t n = t} of [s, t], whose mesh tends 
to zero. Thus, the 1-increment 5 f is the indefinite integral of the 1-increment g. 
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2.2. Young integration. In this section, we will define a generalized integral f f u dg u for 
aC 1 K ([0,T];M nxd )-function /, and a Cj([0, T]; R d ) -function g, with K + 7 > 1, by means of 
the algebraic tools introduced at Section [2TTL To this purpose, we will first assume that / 
and g are smooth functions, in which case the integral of / with respect to g can be defined 
in the Riemann sense, and then we will express this integral in terms of the operator A. 
This will lead to a natural extension of the notion of integral, which coincides with the 
usual Young integral. In the sequel, in order to avoid some cumbersome notations, we 
will sometimes write J st {f dg) instead of f f f u dg u . 

Let us consider then for the moment two smooth functions / and g defined on [0, T]. 
One can write, thanks to some elementary algebraic manipulations, that: 

Jst{fdg)= I f u dg u = f s (Sg) st + / (5f) su dg u = f s (5g) st + J st (5fdg). (11) 

J s J s 

Let us analyze now the term J(5fdg), which is an element of C2(M n ). Invoking Remark 
12.21 it is easily seen that, for s, u,t G [0,T], 

h sut = [S(J(Sfdg))] sut = (Sf) su (6g) ut . 

The increment h is thus an element of C 3 (IR n ) satisfying 5h = (recall that 55 = 0). Let 
us estimate now the regularity of h: if / G Cf([0, T}; R nxd ) and g G CJ([0, T];R d ), from 
the definition it is readily checked that h G Cj +K (R n ). Hence h G ZCj +K (R n ), and if 
k + 7 > 1 (which is the case if / and g are regular), Proposition 12.31 yields that J{5f dg) 
can also be expressed as 

J(5fdg) = A(h)=A(5f5g), 

and thus, plugging this identity into (fTTT) . we get: 

Jstifdg) = f s (5g) st + A st (5f5g) . (12) 

Now we can see that the right hand side of the last equality is rigorously defined whenever 
/ G C 1 K ([0,T];M nxd ), g G Cj([0, T]; R d ), and this is the definition we will use in order to 
extend the notion of integral: 

Theorem 2.5. Let f G Cf([0, T]; R nxd ) and g G C7([0, T}; R d ), with k + 7 > 1. Set 

Jstifdg) = f s (5g) st + A st (6f6g) . (13) 

Then 

(1) Whenever f andg are smooth function, J s t(f dg) coincides with the usual Riemann 
integral. 

(2) The generalized integral J(f dg) satisfies: 

\Jst(fdg)\ < \\f\U\9lM - + c 7 Jf\U\g\U\t - s\^ +K , 

for a constant c 7jK whose exact value is (2 7+K — 1) . 

(3) We have 

n-l 

Jst{fdg) = lim Y)/^ 5g titi+1 , 

where the limit is over any partition H st = {to = s, . . . ,t n = t} of [s,t], whose 
mesh tends to zero. In particular, J s t(f dg) coincides with the Young integral as 
defined in [32] . 
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Proof. The first claim is just what we proved at equation (TT2i) . The second assertion 
follows directly from the definition ffl3l) and the inequality (fTOl) concerning the operator 
A. Finally, our third property is a direct consequence of Corollary 12.41 and the fact that 
3(f b~9) — ~^f^9i which means that 

J(fdg) = \ld-AS\ (f8g). 

□ 

A Fubini type theorem for Young's integral will be needed in the last section of this 
paper. Its proof below is a good example of the importance of Proposition 12.31 and 
Theorem 12.51 



Proposition 2.6. Assume that 7 > A > 1/2. Let f and g be two functions in C7([0, T] : 
R) and h : {(t,s) G [0, T] 2 ; 0<s<t<T}^Ra function such that h(-,t) (resp. h(t, •) ) 
belongs to C A ([t, T];R) (resp. C A ([0, t]; R)J uniformly in t G [0,T], and 

\\h(n, ■) - h{r 2 , OlkicnAn,] < C\n - r 2\ x - (14) 

Then 

pt nr rt r-t 

I I h(r,u)dg u df r = / / h(r,u)df r dg u , < s <t <T. (15) 

J s J s J s Ju 

Proof. Fix s,t G [0,T], with s < t, and divide the proof in several steps. 

Step 1. Here we see that f f r h(r,u)dg u df r is well-defined. Note that we only need to 
show that J h(-,u)dg u belongs to C^([s, T];R) due to Theorem 12.51 
Let ri,r 2 G [s,t], r\ < r 2 , then Theorem 12.51 (2) gives 



h(r 2 ,u)dg u 



n 



h(r u u)dg u 



< 



(h(r 2 ,u) - h(ri,u))dg u 



h(r 2 ,u)dg u 



< lbll 7 (IIM r 2, •) - h{r u •)lloo J [o,n](»"i - $V + c y! x\\h(r 2 , •) - h(r u 0IU,[o,ri](»"i - s) 7+A ) 
+ 11^117 (\\ h ( r 2, 0lloo,[o,r 3 ](''2 - ^i) 7 + c 7 , A ||/i(r 2 , •) || A ,[o,r 2 ] (r 2 - ri) 7+A ) . 



Hence ffT4D implies our claim. The definition of f J h(r,u)df r dg u follows along the same 
lines. 

Step 2. Let Ii st = {t = s, . . . , t n = t} be a partition of the interval [s, t}. Then, according 
to Proposition 12.51 for any v G [0,£) we have 



n-1 



h(t,u)dg u = lim ^^h(t,ti) (6g) 



(16) 



i=0 



Our assumption (|T4j) allows now to take limits in the equation above, so that we obtain, 
for any < s < t < T, 

n-1 

(17) 

i=0 

In order to see that the relation above holds in C 2 ([0, T];R), it is now enough to check 
that both g 1 and q 2 in dHj) are elements of C 2 ([0, T]; R). 



/t n-1 
h(t, u)dg u = lim h(t, U) 5g UM+1 := q 
|n st |— >o 
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However, the fact that q 1 G ^([0, T]; R) can be proved along the same lines as in Step 1. 
The assertion q 2 G ^([0, T];R) can be proved by observing that the limit defining q 2 st do 
not depend on the sequence of partitions under consideration. In particular, consider the 
sequence (vr n ) n of dyadic partitions of [0, T], that is 



7i 



{0 



to < t" < 



< q n = T}, with t\ 



and set, for all s,t G [0,T], ir™ t = n n H (s,t). Then q 2 st = \im r 



for all < s < t < T, and the same kind of arguments as in [6j Theorem 2.2] yield our 
claim q 2 G C 2 A ([0, T]; R). We have thus proved that (Q2D holds in C 2 A ([0, T]; R). 

5^ep 5. From Proposition 12.31 Step 2 and (TT31) we have 

'n-l \ 



/ / h(r,u)dg u df r = lim / /i(r, ti)(g ti+lhr - g ti 

n-l „ t 

= ,Jim 53/ H r iU) (gt l+1 Ar- g tl ) df r 



Ar J 



i=0 
n-l 

lim > 

|n.*|-MJ ^ 



h(r,U)df r (g u+1 - + / h{r,U) (g u+1 Ar - gt i+1 ) df r 



Moreover, thanks to the Holder properties of / and g, we have 

n-l nt i+1 n-l 



£ 

i=0 



- 9ti)dfr 



U 



< C 53^ i+1 — 



,7+A 



i=0 



as in, J — > 0, and thus 



pt rr n—1 / ft \ 

/ / h(r,u)dg u df r = lim Y^l h(r,ti)df r J (# ti+1 - # ti ) . 
./s A |n st |->o . =o \,/ t . y 

Consequently, Step 2 and Theorem 12.51 imply that ([1511 is satisfied and therefore the proof 
is complete. □ 

The following integration by parts and Ito's formulas will be also needed in the last 
part of this paper. 

Proposition 2.7. Let f and g be two functions in Cj([0, T]; R), with 7 > 1/2. Then 



ft9t = fo9o 



/ f u dg u + \ g u dfu, te[0,T]. 
Jo Jo 



Proof. Set q t := f t g t - J f u dg u - / g u df Ul t G [0, T] . It is easy to see that this funcion 
belongs to Cf ([0, T]; R) because of the equalities 

ftgt - f s g s = fs(Sg) st + g s (5f) st + (Sg) st (5f) st 

and 

j-t rt 

f u dg u + / g u dfu = fs(Sg) st + g s (Sf) s t + A st (5fSg) + A st (SgSf), 
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which follows from (fl3|l . Now, since q G C 1 7 ([0, T]; R), with 27 > 1, q is a constant 
function. Otherwise stated, qt = qo = foQo- Therefore the announced result is true. □ 

Proposition 2.8. Let g and h be in C7([0, T],R) and / G C$(R). /et x t = x + 

J^g s dh s , t G [0, T]. T/ien 

/(s t ) = /(so) + / f'(x u )g u dh u , t G [0,T]. 
Jo 

Proof. Proceeding as in the proof of Proposition 12.71 and using the mean value theorem, 
we can show that 

qt = f(x t )- [ f'(x s )g s dh s , te[0,T], 
Jo 

is a 27-Holder continuous function. Therefore the result holds. □ 
Remark 2.9. Proposition 12.81 has been proven in [33] using Riemann sums. 

3. Young delay equation 
Recall first that we wish to consider a differential equation of the form: 

Vt = £0+ / f(Z y u )dx u , t G [0, T], (18) 
Jo 

zi = £■ 

In the previous equation, the integral has to be interpreted in the Young sense of (TT3L the 
initial condition £ is an element of Cj{[— h, 0]; R n ), the driving noise x is in C7([0, T]; R d ), 
with 7 > 1/2. We seek a solution y in the space C| T (IR n ) for 1/2 < A < 7, and / is a 
given function / : C\ ([— h, 0]; R n ) — > IR nxd . In this section, we shall solve equation ([TBI) 
thanks to a contraction argument, and then study its differentiability with respect to the 
driving noise x. Of course, the main application we have in mind is the case where a; is a 
d-dimensional fractional Brownian motion, and this particular case will be considered at 
Section [H 

3.1. Existence and uniqueness of the solution. In order to solve equation (1T81) . some 
smoothness and boundedness assumptions have to be made on our coefficient /. In fact, 
we shall rely on the following hypothesis: 

Hypothesis 1. There exist a positive constant M and X G (1/2,7) such that 

|/(C)| <M, and |/(C 2 )-/(Ci)i<M sup |C 2 (0) - d(0)| 

ee[-h,o] 

uniformly in C, Ci, C2 e C$([-h, 0]; M n )- 

Actually we will assume that / satisfies a stronger Lipschitz type hypothesis on the space 
C^(IR n ). Let us state first a preliminary result before we come to this second assumption: 

Lemma 3.1. Let a = (a l5 a 2 ), with < a\ < a 2 < T, let also Z G Ci([ai — h, a 2 ]; R n ) and 
set 

[Utoz] a = f(Zf), se[ ai ,a 2 ]. 
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Then Hypothesis^ implies that US^> is a map from CiQai — h, 02]; M n ) into ([01,02]; 
M nxd ), satisfying: 

||W (a) Z|L r .<M||Z|k ro hn] . 

II IIA,[ai,a 2 ] — 11 NA,[ai— /i,a2j 

Proof. The proof of this result is an immediate consequence of the definition (JOJ) of Holder's 
norms on C\ and Hypothesis [TJ 

□ 

With this preliminary result in hand, we can now introduce our second hypothesis on 
the coefficient /. 

Hypothesis 2. Taking up the notations of Hypothesis d consider an initial condition 
p G Ci([ai — h, Oi]). We assume that, for any N > 1, there is a positive constant such 
that: 

\\U^{Z x )-U { - a \Z 2 )\\ XM < c N \\Z x - z 2 \\ Uai _ hM , 
for all < Oi < a 2 < T and Z\,Z% G Cp ai a2 {R n ), satisfying 

max { ||^l|U,[ai-/i,a a ]; 11^2 1| X,[ai-h,a 2 ] } < N, 

where A is given in Hypothesis^ 

Observe that Hypothesis [2] holds in particular if, for A > 0, the map admits a 
derivative which is locally bounded, uniformly in a G [0, T]. 

Now that we have stated our main assumptions, the following theorem is the main 
result of this section. 

Theorem 3.2. Under Hypotheses^ and\^ the delay equation [T8\) has a unique solution 
in Cl 0>T (R n ). 

Before giving the proof of this theorem, we establish and auxiliary result. This will be 
helpful in order to get the existence of an invariant ball under the contracting map which 
gives raise to the solution of our equation. 

Lemma 3.3. Let x G C7([ai, a 2 ]; M d ) with 7 > 1/2 and < ai < o 2 , A G (1/2, 7) and v G 
R n . Set a = (01,02), recall notation ffi) and define V (a) : C${[a u a 2 ]; R nxd ) -> C* ai>a2 {R n ) 
by: 

[V (a) Z] s = v + J aiS (Zdx), sG[o!,o 2 ], 
where J aiS (Z dx) stands for the Young integral defined by ^^). Then 

||V (a) Z||A,[ ai ,a 2 ] < ||x|| 7 (||2 , ||oo,[a 1 ,a 2 ]( a 2 ~ Oi) 7 ~ A + C A+7 || Z\\ x,[ai,a 2 ] ( a 2 ~ Ol) 7 ) , 

with c A+7 = (2 a+ t -2)- 1 . 

Proof. Let a x < s < t < T. Then Theorem 12.51 point (3) implies that 

[V {a) Z] t - [V {a) Z] s = J st (Zdx). 

Our claim is then a direct consequence of Theorem 12.51 point (2) and of the definition (J6j) . 

□ 
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Proof of Theorem \3.2c This proof is divided in several steps. 

Step 1: Existence of invariant balls. Let us first consider an interval of the form [0,e], 
which means that, when we include the delay of the equation, we shall consider processes 
defined on [— h, e]. More specifically, let us recall that the spaces C| e (IR ri ) have been 
defined by relation (jHJ). Then we consider a map T : C^ 0e — > C| 0e , where we have set 
C| 0£ = C^ 0e {R n ) for notational sake, defined in the following way: if z G C|o e , then 
T(z) = z, where z\ = & f° r t G [— h,0], and: 

{5z) st = J st {Zdx), with Z u = f{Z z u ), for s,te[0,e]. (19) 

We shall now look for an invariant ball in the space C^ Q e for the map T. 

So let us pick an element z, such that ||z||a,[-a,6] < ^i an d set T(z) = z. On [— h, 0], we 
have z — £, and hence ||<55||a,[-/i,o] = ||^6||a,[-7i,o] = We shall thus choose Aq > 2N%. 

On [0,e], we have now, invoking Lemma 13731 

||^||a,[o, £ ] < ll^l|oo||^|| 7 £ 7_A + c 7iA ||Z|| Ai[0>e ]||a;|| 7 £ 7 . (20) 

Furthermore, according to Hypothesis CD, we have ||Z||oo < M and thanks to Lemma I3TT1 
we also have ||Z|| A ,[o,e] < IN|a,[-m < M N\, by assumption. Then we can recast the 
previous inequality into: 

IN|a,[o j£ ] < M ||x|| 7 e^ A [1 + c^ x N l£ x ] . (21) 

Let us choose now e and Ni in the following manner (notice that e does not depend on 
the initial condition £): 

e = [4Mc 7iA ||x|| 7 ]~ 1/7 A 1, and Aq > 4M||x|| 7 . (22) 

With this choice of e, Aq, inequality ([21]) becomes ||55||A,[o,e] < N\/2. Summarizing the 
considerations above, we have thus found that: 

e = [4Mc 7iA ||x|| 7 ]" 1/7 A 1, Nx > sup {2% 4M||x|| 7 } 

sup {||£2||a,[-m; II^IU,[o,e]} < -y- (23) 

Consider now s < t, with s G [— h, 0] and t G [0, e]. Then, owing to the previous 
relation, we have: 

\(Sz) st \ < \(5z) s0 \ + \(8z) ot \ <^{s x + t x ) <m\t- s\\ 

which, together with the last inequality, proves that .8(0, Aq) in C^ e is left invariant by 
T, under the assumptions of (|23|) . 

Assume now that we have been able to produce a solution to equation ([181) on the 
interval [— h, e]. We try now to iterate the invariant ball argument on [e — h;2e]. The 
arguments above go through with very little changes: we are now working on delayed 
Holder spaces of the form C x {1) e2e , and the map V is defined by T(z) = z, with z = y^- 1 ' 
on [e — h;e\, and Sz having the same expression as in (fill on [e, 2e]. We wish to find a 
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ball .8(0, N 2 ) in C* (1) g 2g , left invariant by the map T. With the same computations as for 
the interval [—h,e], the assumptions of inequality ([23il become: 

e = [4Mc 7iA ||x|| 7 ]~ 1/7 A 1, N 2 > sup {2N yW ] 4M||x|| 7 } . 

Notice again that we are able to choose here the same e as before, by changing Ni into N 2 
according to the value of ||2/^||A,[e-/»,e]- It is now readily checked that B(0, N 2 ) is invariant 
under T, and this calculation is also easily repeated on any interval [ke — h,(k + l)e] for 
any k > 0, until the whole interval [0, T] is covered. 

Step 2: Fixed point argument. We shall suppose here that we have been able to construct 
the unique solution y to (fT8l) on [—h; le], and we shall build the fixed point argument on 
[le — h; (/ + l)e]. On the latter interval, the initial condition of the paths we shall consider 
is = y on [le — h; le]. If T is the map defined on C^ L1 le n + i\ e by fl9jh then we know 
that B(0, iV; + i) is invariant by T. 

In order to settle our fixed point argument, we shall first consider an interval of the 
form [le — h;le + rj\, for a parameter < 77 < e to be determined. On C^i le l£+ri , we define 

a map, called again T, according to (fT9l) . Pick then two functions z Y ,z 2 G C^,i kJe+ , set 

z* = r(z*) for i = 1, 2 and C = z 2 - z 1 . Then ( G Cj te>te+)J , and if le < s < t < le + rj, we 
have 

(SOst = Jst {(Z 2 - Z 1 ) dx) , where Z* = f(Z z *). 
Thus, just like in fl20l) . we have: 

ll^CIUjie-h.fe+jj] < WZ 1 — ^ 2 ||oo,[fe,ie+r ? ]||^|| 7 '7 7 " A + ^^Z 1 — Z 2 \\ Xi [l £j l £+ri ] \\x\\ 7 T) 7 . 

Furthermore, [[Z 1 - ^ 2 ||oo,[k,« £ +»?] < ll^ 1 ~ ^IU.^e+t?] Hence, 

I^ClU.I/e-^ie+r?] < (1 + C 7 ,a) H^ 1 ~ Z 2 \\ X ,[l £ ,le+ V ] \\x\\<y ■ 

We also have Z l - Z 2 = f(Z z ) - f(Z z ), and thanks to Hypothesis El we obtain: 

||*ClU,[fc-h,k-Hj] < (1 + C 7 ,a) ||x|| 7 CAr i+1 ?7 7 H2; 1 - Z 2 \\x,[le-h,le+v}- 

Therefore, we are able to apply the fixed point argument in the usual way as soon as 
(1 + c 7iA )cAT i+1 ||x|| 7 ?7 7 < -, or 77 = [2(1 + c 7iA ) c Nl+1 ||x|| 7 ]~ 1/7 A e. 

With this value of 77, we are thus able to get a unique solution to (fT8l) on [le — h; le + 77]. 

Let us proceed now to the case of [le + 77 — h, le + 2rj\. The arguments are roughly the 
same as in the previous case, but one has to be careful about the change in the initial 
condition. In fact, the initial condition here should be £ 1 ' 2 = y on [le + r\ — h, le + 77] . 
However, we can also choose to extend this initial condition backward, and set it as £ />2 = y 
on [le — h, le + 77]. We then define the usual map T as in (1T9I) . and we have to prove that 
B(Q,Ni + i) is left invariant by T. To this purpose, take z G C^, 2 l£+v le+2r) in £>(0,iV; +1 ), 
and set z = T(z). Observe then that, for any t G [le + 77, le + 277], we have 

h = Cl +V + f fiK) dx u = 4+ f £+V f{Zl) dx u + f f(Z z u ) dx u = Cl+ f f{Z z u ) dx u , 

Jle+rj Jle J le+rj J le 

where we have used the fact that £ 1 ' 2 = y on [le — h,le + 77] solves (ITSil . It is now easily 
seen that z is in B(0, iVj+i), and this allows to settle our fixed point argument as in the 
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previous case, with the same interval length rj. This step can now be iterated until the 
whole interval [le; (I + l)e] is covered. 

□ 



3.2. Moments of the solution. The moments of the solution to (fT8l) can be bounded 
in the following way: 

Proposition 3.4. Under the same assumptions as in Theorem \3.S\ let y be the solution 
of equation ( fiffj] on the interval [0, T], with an initial condition!; G C A ([— h, 0]; IR n ). Then 
there exists a strictly positive constant c = 0(7, A, M, T) such that 

\\y\\x,[-h,T\ < cmax [U\\x, M^^, \\x\\ 7 ] . 

Proof. From the proof of Theorem 13.21 we know that \\y\\ \,[-h,T\ is finite. Let us assume 
that this quantity is equal to A, and let us find an estimate on A. One can begin with 
a small interval, which will be called again [0, e], though it won't be the same interval as 
in the proof of Theorem I3.2L In any case, taking into account that y solves equation (TT8]) . 
we obtain similarly to (1201) : 

INU,[o, e ] < M\\x\\ y e^ x + c y , x M\\5y\\ x , Hh£] \\x\\^ 

< M ||a;|| 7 e 7 - A + c 7jA MA ||x|| 7 e 7 = g(e,K). (24) 
Along the same line, for any k < [T/e], we have 

||^||A,[te,(fc+l)e] < g(£, K )- 

Take now s,t G [0, T] such that ie < s < (i + l)e < je < t < (J + l)e. Set also = s, 
tfe = ke for z + 1 < k < j, and tj + x = t. Then 



\(Sy) 



.sf 



< A') - 4) A < <7(e, A)(j - i + l)^ A (t - s) 



k=i 



where we have used the fact that r \— > r x is a concave function. Note that the indices i, j 
above satisfy (j — i + 1) < 2T/e. Plugging this into the last series of inequalities, we end 
up with 



¥v\\w 



T 



< 



9(e,K)(2T) 

c-l-A 



1-A 



M ||x|| 



c^ x M K\\x\Ue> +x - 1 



(2T) 



l-A 



Thus the parameters A and e satisfy the relation: 

M \\x\ 



A < 



n 



In order to solve (1251). choose £ such that 



c J)X MK\\xh£' r+x - 1 



c 7iA M||a;|| 7 e 7+A - 1 (2T) 1 - A 



(2T) 



l-A 



lieii; 



(25) 



that is 



e= [2c 7 ,aM||x|| 7 (2T) 1 - a ]" 1/(7+A ~ 1) . 

Plugging this relation into ([251) . we obtain the result when e <T. 

Finally, T < £ if and only if T 7 < [2 2 ~ A c 7+A M| |x| ^j" 1 . Thus, by inequality (El), the 
proof is complete. 

□ 
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3.3. Case of a weighted delay. In this subsection, we prove that our Hypotheses [Tj 
and [2] are satisfied for the weighted delay alluded to in the introduction, that is for the 
function / given by equation (J2J). 

Proposition 3.5. Let v be a finite measure on [— h,0] and a : W 1 —> W nxd a four times 
differentiable bounded function with bounded derivatives. Then Hypotheses U\ and [H are 
fulfilled for f : C?([-h, 0]; R n ) -> R nxd defined by: 

f{Z)=o(J Z{Q)v{d6) 

with Z e C£([-h, 0];R n ). 

Proof. We first show that Hypothesis Q] holds. More specifically, the condition |/(C)| < M 
being obvious in our case, we focus on the second condition of Hypothesis [TJ Let Z\,Z 2 G 
C^([-h, 0};R n ). Then there is a constant C > such that 

1/(^0-/(^)1 
-o 



< C 



[Z 1 (9)-Z 2 (0))u(d9) 



<Cv([-h,0]) [ sup \Z 1 (6)-Z 2 

6e[-h,0] 



Therefore Hypothesis [U is satisfied in this case. 

Now we prove that is Frechet differentiable in order to analyze Hypothesis El Since 
the map Z t— > J,Z(- + 9)v(d6) is easily shown to be a bounded linear operator from 
C£([ai - h,a 2 ];R n ) into C^([a 1 ,a 2 \;R n ), we only need to show that 

o-:C^ a2 (R n )^Cl ai! jR n * d ), with p±a(p), 

is Frechet differentiable in the directions of Co aia2 {R n ), with derivative [Dcr(Z)£](t) = 
a'(Z(t))£(t). Towards this end, we have to show that, taking Z G Cp ai a2 {R n ) and £ G 
C A ai ,a 2 (K n ),and setting 

q t = a(Z(t) + t{t)) - a(Z(t)) - a'(Z(t)) £(t), 

then 

hm jfe^ = 0. (26) 

Uh,[ ai -h,a 2 ]->0 FlU,[ai-/i,a 2 ] 

In order to prove relation (l26l) . define a function b : [0, l] 2 — > R by: 

6(A, fj) = Z(s) + X£(s) + n[Z(t) - Z{s)\ + A// [£(t) - £(s)} . 

Observe then that 6(1, 1) = Z(t) + £{t), 6(1, 0) = Z(s) + £(s), 6(0, 1) = Z(t) and 6(0, 0) = 
Z(s). We will also set H(\,fi) = o~(b(\, //)). Then 

a(Z(t) + £{t)) - a(Z(t)) - a'{Z{t)) £{t) 
= (7(6(1, 1)) - (7(6(0, 1)) - (7'(6(0, 1))[6(1, 1) - 6(0, 1)] 
and similarly, we have: 

a{Z{s) + £(s)) - o-(Z(s)) - a'{Z{s)) £{s) 



l -j\l x H{\l)[l-X]d\ 



[ dl x H(\,0)[l-\] dX. 
Jo 
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Hence, plugging these two relations in the definition of q, we end up with: 

(5q) st = [ (o^ x H(X, 1) — d1 x H(X, 0)) [1 — A] dX 
Jo 

= [ d? XXfl H(X, 0) [1 - A] dX + [ di x ^H(X, /*) [1 — A][l — fi] dXdf,. 
Jo J[o,ip 

The calculation of d xx H(X,0) and d^ x H(X,/i) is a matter of long and tedious compu- 
tations, which are left to the reader. Let us just mention that both expressions can be 
written as a sum of terms from which a typical example is: 

a"'(b(X, fi)) l(5Z) st + fi(5Z) st ] [£{s) + X{6£) st ] (6£) st . (27) 

These terms are obviously quadratic in £, and can be bounded uniformly in A, fi, s, t 
under the hypothesis a G C%. Notice that, in order to bound the term \£(s)\ in ([27]) . we 
use the fact that £ has a null initial condition, which means in particular that \£(s)\ < 
(02 — cti + A H^ll a,[q,i— /i,a 2 ] • This finishes the proof of (1261) . The continuity of Dcr(Z) and 
the existence of the constant introduced in Hypothesis [2] are now a question of trivial 
considerations, and this ends the proof of our proposition. 

□ 

Remark 3.6. The proof of Frechet differentiability of / was not necessary for the existence- 
uniqueness result, which relied on some Lipschitz type condition. However, this stronger 
result turns out to be useful for the Malliavin calculus part, and this is why we prove it 
here. Nevertheless, notice that Theorem 13.21 holds true for a C% coefficient a. 

3.4. Differentiability of the solution. In this section we study the differentiability of 
the solution of ([TBI as a function of the integrator x, following closely the methodology of 
|26| . In particular, our differentiability result will be achieved with the help of the map 

F : Cj AT (M d ) x C^ Q<T (R n ) - C^ T (R n ) given by 

[F(k, Z)] t = Z t - Jot (f(Z Z+i ) d(x + fc)) , te [0, T] (28) 

where £t = £o for t G [0, T], and £t = £t for t G [—h, 0]. Here we recall that £ stands for an 
initial condition in C±([—h, 0]). In this section the coefficient / will satisfies the following: 

Hypothesis 3. Set t = (0,t), and recall that the map U^' has been defined at Lemma 
[Q We assume thatU^ : C^W 1 ) -> C A ([0, t\\ R nxd ) is continuously Frechet differen- 
tiable in the directions ofC$ At (R n ), for some X G (1/2,7). We call VW (t) : C| At (R n ) -> 
£(C^ t {R n )]C^ 0>t (R nxd )) its differential, where C{C^ t (R n );C^ t (R nxd )) denotes the lin- 
ear operators from C^ o t (R n ) into C^ 0tt (R nxd ). Moreover, for s < t and Z G C^ 0T (R n ), 

\VU^\y)]{Z) = [VU^\y)}{Z) on [0,s], 
where y is the solution of equation (T8\). 

Remarks 3.7. (1) Notice that we have shown, during the proof of Proposition 13.51 that 
the weighted delay given by ((2l) also satisfies this last assumption. 

(2) If Z G C A o t (M rt ), then 

||W«(i/)(Z)|| A , M <|VZ^(y)|||Z|| A ,| Pl i|. 
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Indeed, set Z s = Z s for s G [0,£], and Z s = Z t for s > t. Therefore Hypothesis [3] implies 
\\VUV(y)(Z)\\ xm < \\VU^(y)(Z)\\ mT] < \VU^\y)\\\Z\\ um =< \VU^(y)\ \\Z\\ xm , 
and our claim is satisfied. 



We are now ready to prove the differentiability properties for equation (ITBj) : 

Lemma 3.8. Under the Hypothesis^ the map F given by $28\) is continuously Frechet 
differentiable. 

Proof. Let us call respectively D\ and D 2 the two directional derivatives. We first observe 
that, for k,g G Q AT {R d ) and Z G C^ 0T {R n ), we have: 

F(k + g,Z)-F(k,Z) + [ \u^\Z + i)] dg s = 0. 

Jo 1 Js 

In other words, the partial derivative D\F is defined by 

D 1 F(k, Z)(g) = - J \l4^\Z + 0] s dg s = -J . {[U (T \Z + £)] dg) . 

We shall prove now that DiF is continuous: consider k,k G CQ 0T {R d ) and Z,Z G 
CQ 0T (M ri ). For notational sake, set also || ■ ||a for || ■ ||a,[o,t]- Then, according to Lemma lcTBl 
we obtain: 



D 1 F{k,Z){ri)-D x F{k,2!){r{) 

<Hl t {\u {T) {z + l)-u^\z + l) 



j([u( T \z + o-uV\z + 0}d V( 

rp-y-X 



u^(z + i)-u^(z + i) 



which, owing to Hypothesis [3J implies that DiF is continuous. 

Concerning D 2 F we have, for k G Cj 0T (M d ), Z G C A 0T (M n ) and Z G C A 0T (M n ), and 
thanks to Theorem 12.51 



F(k, Z + Z)- F(k, Z)-Z + J {[VU {T \Z + 1)]{Z) d(x + fc)) 

fe|| 7 (||w (T) (z + z + £)- u {T \z + e) - [vw (T) (z + £)](£) 



< \\x 



T 



7-A 



+ z + o - w (T) (^ + - [w (T) (z + £)](£) 

Therefore, making use of Hypothesis 02, we have that: 

D 2 F(k,Z)(Z) = Z- [ VU^(Z + 0(Z) s d(x s + k s ). 

Jo 

The continuity of D 2 F can now be proven along the same lines as for DiF, and the 
computational details are left to the reader for sake of conciseness. The proof is now 
finished. 

□ 



The following will be used to show that D 2 F(k, Z) is a linear homeomorphism. 
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Lemma 3.9. Let w G Cq 0T (M. u ), y the solution of flfy) and assume Hypotheses^, and 
hold. Then the equation 

Z t = w t + [ ([VU^ T \y)}(Z)) s dx s , 0<t<T, (29) 
Jo 

has a unique solution Z in Coor(^ n )- 



Proof. Similarly to the proof of Theorem 13 .21 we choose e G (0, T) and set T : Cq 



A /Dill 



15771(2) yield 



given by T (Z) = w + J .([W(^'(y)}(Z) dx). Then, Lemma 1531 and Remark 



%{Z) - %{Z) 



A,[0,e] 



J ([VU {£ \y)}(Z - Z) dx 



A,[0,e] 



7-A 



VU^{y){Z-Z) 



oo,[0,e 



+ c A+7 T ; 



vu&(y){z-z) 



A,[0,e] 



< |VW (T) (y)k 7 - A ||x|| A ||Z - Z||U,[o, £ ](r A + c A+7 T A ). 
That is, for e small enough there exists < C < 1 such that 



%{Z)-%{Z) 



A,[0,e] 



< c\\z-z\ 



|A,[0,e] 



Hence, by standard contraction arguments, one can find a unique Z e G Cq (K n ) such 
that 

^ £ = w t + [ ([VU^(y)}(Z s )) s dx s , 0<t<e. 
Jo 

Now we introduce % : C^^ 2e (R n ) -> C| £ e 2e (M n ) defined by 

f e (Z)(t) = w t -w e + ZI + J ([VU^\y)}(Z)) s dx s , t G [e,2e]. 

Then, as in the beginning of this proof, we have 

%{Z)-%(Z) <C\\Z-Z\\\ HQM . 

A, [0,2£j 

Therefore, there is a unique 

Z 2e e C^ e e 2e (M n ) such that 



f \[VU^\y)]{Z^)) s dx sl 0<t<2e, 
Jo 



due to Hypothesis [3j 



Finally by induction, we can figure out a function Z ke G C^ (k _ 1)e ( fc _ 1 ) e j fce (^ n ) sucri that 



[ {[VU^\y)](Z ks )) s dx s , 0<t<ke. 
Jo 



Consequently, by Remark l3"771 (2). it is not difficult to see that Z t = Z\ e for t G [(k— l)e, fee] 
is the unique solution to equation (1291) . □ 
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Proposition 3.10. Assume that Hypotheses^ to\E are satisfied. Let y be the solution of 
equation ( TTgj) . Then the map h i— > y(x + h) is Frechet differentiate in the directions of 
Cj T (]R d ) ; as a C| T (]R n ) -valued function. Moreover, for h,k G Cj or (IR d ) ; we have 

[Dy(x)(k)] t = fu^\y{x)) s dk s 
Jo 

+ / [VU^\y(x))(Dy(x)(k))} s dx s . (30) 
Jo 

In particular, [Dy(x)](k) is an element ofCQ 0T (M. n ). 

Remark 3.11. Let us recall that equation f|30l) has a unique solution, thanks to Lemma lcT9l 

Proof of Propositio ns. 1 (k Like in [26], the proof of this result is a consequence of the 
implicit function theorem, and we only need to show that D 2 F(0,y(x) — £) is a linear 
homeomorphism from CQ 0T (IR n ) onto CQ 0T (IR n ). Indeed, in this case we deduce that 
h i— > y(x) is Frechet differentiable with 

Dy{x)(k) = - (D 2 F{h,y{x) - I))' 1 o D x F{h,y{x) -£)(k), (31) 

which yields that (l30l) holds. 

Finally, notice that D 2 F(0,y(x) — £) is bijective and continuous according to Lem- 
mas 13.81 and 13.91 Consequently the open mapping theorem implies that the application 
D 2 F(0, y(x) — £) is also a homeomorphism. 

□ 

Interestingly enough, in the particular case of the weighted delay of Section 13.31 one 
can also derive a linear equation for the derivative [Dy(x)] t , seen as a Holder-continuous 
function. 

Proposition 3.12. Let a and u be as in Proposition \ 3.b~[ Let also f and y be defined 
by (HP and (Q2P. respectively. Assume that v is absolutely continuous with respect to the 
Lebesgue measure with Radon-Nykodim derivative in L p ([—h,0]) forp > 1/(1 — 7). Then, 
for i G {1, . . . , n} and k G Cq 0T {R u ), we have 

Dyl(x)(k)=J2 f^{r)dkU 
3=1 Jo 

where, for j G {i, . . . , d} and i G {1, . . . , n}, $ tJ ' is defined by the equation 

n d „ t 

*?(r) = (U^\y)y> + {([VU^(y)ry l (^(s))) s dx l s , < r < t < T, (32) 

m=l 1=1 ^ r 

and $ t (r) = for all < t < r < T. 

Remark 3.13. Note that, for each s G [0, T] equation (132]) has a unique solution in 
C A ([s,T];M n ) due to Lemma ITT91 
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Proof of Proposition ^. 12l In order to avoid cumbersome matrix notations, we shall prove 
this result for n = d — 1: notice that an easy consequence of the proof of Proposition 13.51 
is that in our particular case, 

[vu^\z){k)] t = a' (J° z^eu(de) \ (J^kt+guide)) . (33) 

Set now q t = o~(f\yt+o v(d6)) and q' t = a'(f\y t+ gi'(d8)), and write y = y(x). Then 
equation fl30l) can be read as: 

[Dy{k)] t = [ q s dk s + U t , with U t = f q' s ( [ [Dy(k)] s+e u(d9)) dx s . (34) 

Jo J0 \J-h J 

The Fubini type relation given at Lemma I2U1 allows then to show, as in [26], Proposition 
4], that 

[Dy(k)] t = / $ t (r)dkr, (35) 
Jo 

for a certain function <E>, A-Holder continuous in all its variables. In order to identify the 
process $, plug relation (1351) into equation (1341 and apply Fubini's theorem, which yields 



U t = I v{d6) I q'A I $ s+e {r)dk r I dx s 



It should be noticed that this point is where we use the fact that v(d9) = fi(9) d6 with 
G L A ([— r, 0]). Indeed, in order to apply Lemma [2T61 to x, k and r] i— > F{rj) = f_, fi(9) d9, 
we will assume (though this is not completely optimal) that F is 7-Holder continuous. 
However, a simple application of Holder's inequality yields 

\F(r)2) - F( Vl )\ < c\t - s\^ p |M|lp([-M)- 

It is now easily seen that the condition {p — l)/p > 7 imposes p > 1/(1 — 7). 

Owing now to a (slight extension of) Lemma \2. 61 we can write 

ro r(t+8)+ rt 

U t = v{dQ) I m t (r,9)dk r , with m t (r,9) = / q' s § s+e (r) dx s . 

J-h Jo Jr-0 

Apply Fubini's theorem again in order to integrate with respect to k in the last place: we 
obtain 

U t = f ([ m t {r,9)v(d9)) dk r = [ ([ u{dd) [ q' s <f> s+e (r) dx s ) dk r , 

Jo \J-[(t-r)Ah] J Jo \J-[(t-r)Ah] Jr~9 J 

and going back to fl34l) . which is valid for any A-Holder continuous function k, we get that 
$t is defined on [0, t] by the equation 

$t(r) =qt+ [ ( I q' s $ s +e(r) rfx s ) v(d8), 

J-[(t-r)Ah] \Jr-e J 

and $t(r) = if r > t. A last application of Fubini's theorem allows then us to recast the 
above equation as 

$t(r) = q t + / q' s ( [ $ s +e(r) u{dB)) dx s . 

Jr \J-[hA(s-r)] J 
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Notice now that, if 6 < — (s — r) in the above equation, then s + 6 < r, which means that 
$ s+ e(r) = 0. Hence, we end up with an equation of the form 



Mr) = Qt + q'sij ®s+e( r ) v(d6) \ dx s 



which is easily seen to be of the form (|32l) . 

□ 

3.5. Moments of linear equations. In order to obtain the regularity of the density for 
equation (fT8l) . we should bound the moments of the solution to equation (1291) . This is 
obtained in the following proposition: 

Proposition 3.14. Let f be a mapping from C^ 0T (R n ) into the linear operators from 
C A 0T (M n ) into C x {[0,T];W nxd ) such that, for < a < b < T, y e C 5 A or (M n ) and z G 

oo,[a,b] oo,[a— h,b] ■ 

(2) \\f(y)z\\x,[a,b] < M\\z\\x,[a-h,b] + M\\y\\x,[a-h,b]\\z\\oo,[a-h,b]- 

Also let y be the solution of the equation (Q2P, w G CQ 0T (IR n ) and z G CQ 0T (IR n ) the 
solution of the equation 

z t =w t + [ (f(y)z)(t)dx t , te[0,T\. 



Then 

\\4x,[o,t] < Ci||w|| Ai[ o,r]^, A e C2D7 ' A , 
for two strictly positive constants q = Cj(T, 7, A, M), i = 1, 2 and 

^ A = (UhWxh) 1 ^ + \\x\\yi + || x ||PA+7-l)/((7+A)(7+A-l))_ 

Remarks 3.15. (1) Observe that if / is as in Proposition 13.51 and / = WU^ T \ then 
straightforward calculations show that Conditions (1) and (2) in the Proposition are 
satisfied. 

(2) The fact that zo = implies that 

INU,[o,T] < cxTiHIw]^^ 2 ^. 

(3) Let A = 7. Then (7 + 2A - l)/((7 + A) (7 + A - 1)) in Proposition I3T41 is smaller than 
2 for 7 > H , where H = (7+ vT7)/16 ~ 0.6951. This is the threshold above which our 
general delay equation will admit a smooth density. 

(4) The unusual threshold H above stems from the continuous dependence of the solu- 
tion on its past, represented by the measure v. In case of a discrete delay of the form 
a(y t ,y t _ ri , . . . ,y t - r ), we shall see that all our considerations are valid for any H > 1/2. 



Proof of Proposition [X7^[ We first consider two generic positive numbers k G N and e, 



such that (k + l)e < T. Then Theorem 12.51 point (2), and Conditions (1) and (2) imply 

\\ z — HU,[fc£,(fc+l)e] 



< \\f(y)z\\ OC:[k£:{k+1)£] \\x\\ J e' y ~ + c Jt x\\f(y)z\\x,[ke,(k+i)e]\\x 

< M\\z\\oo,[0,(k+l)e]\\ x \U 
+C 7jA M||a;|| 7 (||z||A,[0,(fc+l)e] + ||2||oo,[0,(*+l)e]||j/|U,[0,T|) e" 



7-A 
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The following (arguably non optimal) bound on ||^||oo,[o,(fe+i)e] can now be easily verified 
by induction: 

fc+i fe+i 

|H|oo,[0,(fc+l)e] < yj2 fc+1 ~ t ||z — Z(t-l) E ||oo,[(i-l)e,ie] - 2 k+1 ~ % \\ z\\ \\{i-l) e ,ie\ ■ 

i=l i=l 

This yields 

\\Z - W\\x\ke,{k+l)e\ 

/k+1 

< M\\x\\^ I ^2 k+1 ~%z - z {t _ 1)£ \\ U[l _ l)eM 
+c 7>A M||x|| 7 £ 7 (H^IUjo.te] + IMU,[fc £ ,(/c+i) e ]) 

+c 7 , A M||x|| 7 ||y|| Ai[0iT] ^ +A ^2 fc+1 -i^-^_ 1)£ || A , [(l _ 1)£ , £] j . (36) 

Now the proof can be split in three steps. 
Step 1. Bounds depending on e. Let 

e = (T + [6Af||x|| 7 (l + c 7 , A )]^ + [eMUxll.c^lblU,^]] 17 ^)" 1 A T. (37) 
Note that in this case, inequality (1361) yields 

\\z\\\,{ke,(k+l)e] 

k 



, i=l 

k 



+c 7iA M||x|| 7 ^ 2||z|| A|[0l te] + e%\\ XM J2 2k+2 ~l z h^-i)e,ie] 



i=l 



< 2\\w 



X,[ke,{k+l)e] 
k 

+ Y^2 k+2 ~%z\\ m . l)eM (M||x|| 7 e^ + ^ )A M||x|| 7 e^ + c 7 , A M||x|| 7 e'^ A || 1 /|| Ai[0> r|) 



i=l 



< 2|Mk[fce,(fc+i) e ] + ^2 fc+1 i |k|| A ,[(j-i) £ ,j £ ], (38) 

i=l 

where we have used (l37l) in the last step. 

Step 2. Bounds for ||^|| A ,[fc e ,(fe+i) e ]. Here we will use induction on k to show that 



MUa-DeM < E 22i+1 ~ 2i IHk [( ;-i)^]- ( 39 ) 

By (1381) we have that this inequality holds for i = 1. Therefore we can assume that 
fl39l) holds for any positive integer i less o equal than k to show that it is also true for 
i = k + l. 
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The inequalities (1381) and (1391) lead us to write 

\\ z \\x,[ke,(k+l)e] 

k i 

< 2 1| ^ || ^^^^fc^^^^] S^^ 1 — ^ S 3 ^ 1 "^ || ^7 1| ^^^^i^^^^] 

i=l 3=1 

k k 

< 2||w|| Ai[tei(fc+1)6] + £ \\w\\ x>[(j - l)e j E] 2 k+2 - 2 i 21 



k 

< 2\\w\\x,[k e ,(k+l)e] + ^2 IHklO-lKje] 



i=l 



92fc+3-2j 



Now it is easy to see that fl39l) also holds for i — k + 1. 

.S^ep 5. Final bound. Let fco such that fco£ < T < (A; + Then, by Step 2 we have 

ikiiw^ikiu^EE 22 ^' 

fc=l j=l 

< \\w\\ mT] (k ) 2 2 2k ^ < \\w\\ um (2T/e) 2 2 2T ^\ 

Thus the proof is finished by plugging relation (1371) into the last expression, and invoking 

□ 



Proposition I3.4L 

The following result is a slight extension of Proposition 13.141 allowing to take into 
account the case of constant but non vanishing functions. 

Corollary 3.16. Let f , -D 7 ,a ; w and y be as in Proposition \3.14\ Furthermore, assume 
that f is a mapping from C^ T (K n ) into the linear operators from the constant functions 
on [-h,T] into C A ([0, T}; W> x ' d ) satisfying the Conditions (1) and (2) of Proposition\3J^ 
when z is a constant function. Then the solution of the equation 



z t = c + w t + I (f(y)z)(t)dx t , t e [0, T], 

satisfies the inequality 

IMU,[o,r] < ci 



w + 


I 







X,[0,T] 



D 2 , e C2D ->' x 

7,A ° : 



where c stands for the constant function c t = c. 



Proof. The proof is an immediate consequence of Proposition 13.141 Indeed, we only need 
to observe that 



z t -d t = w t + / (f(y)c)(t)dx t + / (f(y)(z - c))(t)dx t , t e [0,T], 
Jo Jo 

where c(i) = c, t G [0,T]. □ 
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4. Delay equations driven by a fractional Brownian motion 
Here we consider the Young stochastic delay equation 

Vt = £o + J f{Z v t )dB u 0<t<T, 

Zl = e, (40) 

where B = {B t ; < t < T} is a <i-dimensional fractional Brownian motion (fBm) with 
parameter H G (1/2,1). The coefficient / satisfies Hypotheses [TJ13] and £ is a given 
deterministic function in Cj([— h, 0]; IR n ), for some A < 7 < H . Remember that A G 
(1/2, H) is introduced at the beginning of Section[3l 

The fBm B is a centered Gaussian process with the covariance 

R H (t, s)5 id = E{B\Bi) = \5 id {s 2H + t 2H - \t - s\ 2H ). 

In particular, B has //-Holder continuous paths for any exponent v < H. Consequently, 
from Theorem 13.21 and Hypothesis [TJ131 equation ([401) has a unique C^ 0T (M n )-pathwise 
solution. 

Here, our main goal is to analyze the existence of a smooth density of the solution 
of equation (1401 . This will be done via the Malliavin calculus or stochastic calculus of 
variations. 

4.1. Preliminaries on Malliavin calculus. In this subsection we introduce the frame- 
work and the results that we use in the remaining of this paper. Namely, we give some 
tools of the Malliavin calculus for fractional Brownian motion. Towards this end, we 
suppose that the reader is familiar with the basic facts of stochastic analysis for Gaussian 
processes as presented, for example, in Nualart [23] . 

Henceforth, we will consider the abstract Wiener space introduced in Nualart and 
Saussereau [26], in order to take advantage of the relation between the Frechet derivatives 
of the solution to equation ([4111) (see Proposition 13. 101) and its derivatives in the Malliavin 
calculus sense (see [23], Proposition 4.1.3). This abstract Wiener space is constructed as 
follows (for a more detailed exposition of it, the reader can consult [26]). 

We assume that the underlying probability space (Q, P) is such that Q is the Banach 
space of all the continuous funtions C ([0, T]; M. d ), which are zero at time 0, endowed with 
the supremum norm. P is the only probability measure such that the canonical process 
{Bt] < t < T} is a d-dimensional fBm with parameter H G (1/2, 1) and the cr-algebra 
T is the completion of the Borel cr-algebra of Q with respect to P. 

Two important tools related to the fBm B are the completion H of the M. d - valued step 
funcions S with respect to the inner product ((lp,^], • ■ • , l[o,t d ]), (l[o,si], • • • , l[o,s d ])) = 
J2f=i B-H{si,U) and the isometry K* H : H — > L 2 ([0, T] d ), which satisfies 

^((l-io.ti]) • • • , Mo,t d ]) = (l[o,ti](-)^fl"(*i; •)>•••> l[0M K H(td, ■)), 
where Kn(t, s) = ChS 1 ^ 2 ~ h f (u — s) H ~ 3 ^ 2 u H ~ 1 ^ 2 du is a kernel verifying 

ptAs 

R H (t,s) = / K H {t,r)K H {s,r)dr. 
Jo 
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It should be noticed at this point that K* H can be represented in the two following ways: 
[K* B if]t = j\ r d r K(r,t)dr = ^^2-^-1/2^-1/2^^ (41) 

where Z£L stands for the fractional integration of order a on [0, T] (see [24] for further 
details). 

The isometry K* H allows us to introduce the version of the Reproducing Kernel Hilbert 
space Tin associated with the process B. Namely, Let JCh be given by 

K H : L 2 ([0,T];R d ) ^ H H := IC H (L 2 ([0,T];R d )), (JC H h)(t) = [ K H (t, s)h(s)ds. 

Jo 

The space TC is continuously and densely embedded in fl. Indeed, it is not difficult to see 
that the operator IZr '■ TC — > Tin defined by 

n H <t>= [ K H (; S )(K* H <f))(s)ds 
Jo 

embeds TC continuously and densely into fl, because, as it was pointed out in [26], 1Zh{4>) 
is //-Holder continuous. Thus, we have that (fl, TC, P) is an abstract Wiener space. 

Now we introduce the derivative in the Malliavin calculus sense of a random variable. 
We say that a random variable F is a smooth functional in S if it has the form 

F = f(B(h 1 ),...,B(h n )), 

where hi, . . . ,h n G TC and / and all its partial derivatives have polynomial growth. The 
derivative of this smooth fuctional is the H- valued random variable given by 

VF = Y^^-{B{h 1 ),...,B{h n ))h l . 

For p > 1, the operator V is closable from L p (fl) into L P (Q; H) (see [23]). The closure of 
this operator is also denoted by V and its domain by D 1,p , which is the completion of S 
with respect to the norm 

\\F\\l P = E{\F\ P ) + E(\\VFr H ). 

The operator V has the local property (i.e., VF = on A C fl if IaF = 0). This allows 
us to extend the domain of the operator V as follows. We say that F G if there is 
a sequence {{fl n , F n ),n > 1} C T x D 1 ^ such that fl n f fl w.p.l and F = F n on fl n . In 
this case, we define VF = VF n on fl n . 

It is known that, in the abstract Wiener space (fl,H,P), we can consider the dif- 
ferentiability of random variable F in the directions of H. That is, we say that F is 
7i-difTerentiable if for almost all uj G fl and h G H, the map e ^ F{uj + elZHh) is differen- 
tiable. The following result due to Kusuoka [H] (see also [23], Proposition 4.1.3) will be 
fundamental in the study of the existence of smooth densities of the solution of equation 

(SOD. 

Proposition 4.1. Let F be an Ti-differentiable random variable. Then F belongs to the 
space D^, for any p > 1 . 
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We will apply this result to the solution of equation (T4"0l as follows. Note that for 
if G Tt, we have the inequality 

\(n HV y(t) - (K H tpy(s)\ = (e[\b\ - bi\ 2 }) 1/2 y\\ n < y\\ H \t - S \ H . 

Consequently, Proposition 13.101 (see also Lemma 14.21 below) implies that the random 
variable y t defined in equation (1401 is also 7Y-differentiable, which, together with Propo- 
sition HIH yields that y\ belongs to for every t G [0, T], p > 1 and i G {1, . . . ,n}. 
Moreover, the relation between the 7i-derivative and V is given by (see also Lemma l4~3l) . 

(Vyl h) H = ±-y\{u> + eH H h) | e=0 , h G H. (42) 



More generally, if oj i — > X(oS) is infinetely Frechet diferentiable in the directions of 
Co OT (R), then for a smooth random variable X, then 

(V n X, h x ® • • • ® 

d d 

= D-R,Hhx,...,K H h n X = t^- . . . -^X(uj + eiTl hl + . . . + e n TZ hn )\ £l= „- £n= o. 

4.2. Existence of the density of the solution. In this section we establish that, for 
each t G [0, T], the random variable y t introduced in equation (1401 has a density. 

Let us start with two important technical tools. The first one relates the derivative of 
the vector-valued quantity y t with the derivative of y as a function. 

Lemma 4.2. Let y be the solution of (JU\) and t G [0,T]. Then almost surely, h \— > 
y t (B + h) is Frechet differentiable from CQ 0T {R d ) into MJ 1 . Furthermore 

Dy t (B)(h) = [Dy(B)(h)] t . 
Proof. The proof is an immediate consequence of 

\y t (x + h)-y t (x)-(Dy(x)(h))(t)\ 
= \y t (x + h) - y t (x) - (Dy(x)(h)) (t) 

-y {x + h)- y (x) - (Dy{x){h)) (0)| 
< \\y(x + h)-y(x)-Dy(x)(h)\\ x t\ 

with x,he C$ 0T {R d ). 

□ 

Lemma 4.3. Let y be the solution of jj^j . Then y\ belongs to for every t G [0,T] 
and i G {1, . . . , n}. Moreover, for h EH, we have 

(Vyih) H = [Dy\B){K H h)] t . (43) 

Proof. By Proposition 14.11 and Lemma 14.21 we have already shown that y\ is in for 
every t G [0, T] and % G {1, . . . , n}. 

Furthermore, by (|42l) and Lemma 4.2, we have 

(Vyl h) H = Dn Hh yl = Dy\{B){n H h) = {Dy\B)(K H h)) (t). 

Thus, the proof is complete. □ 

We now use the ideas of Nualart and Saussereau [26] to state one of the main results 
of this section: 
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Theorem 4.4. Let us assume that HypotheseslMEhold, recall that£ is the (functional) ini- 
tial condition of equation (JU\ ), and assume that the space spanned by {(/(O^ • • • > /(£) nj ')i 
1 < j < d} is W 1 . Then for t G (0, T], the random variable y t given by [Jty j is absolutely 
continuous with respect to the Lebesgue measure on IR n . 

Proof. As in [26] (proof of Theorem 8), we have that y\ belongs to D^. Therefore we 
only need to see that the Malliavin covariance matrix 

Q? := (Vyi,Vyi) n (44) 

is invertible almost surely. 

For »el", following [26] (proof of Theorem 8), we have 



T 

v 



Q t v = \(Dy(B)(n H h m )(t),v) Rn \ 



m=l 



where {h n ,m > 1} is a complete orthonormal system of TL. 

Now assume that the Malliavin matrix Q t is not almost surely invertible. Then, on the 
set of strictly positive probability where Q t is not invertible, there exists u 6 M™, d / 
such that v T Q t v = 0. Moreover, recalling our notation ([281) . it is clear from equation (I3B 
that D 2 F(k, Z) is a linear homomorphism. Hence, we obtain that 

= (DiF(0,y(S-|))(7e H /i ro )(t),t;o)R« 

U^ T \y(B)) s dn H h m (s),v 

n d 







EE l, o (U^{y{B))f s dn H y m {s) 

i=l j=l J ° 

n 

(U<V{y(B))y l [0 ,t],h m ) H , for all m > 0, 



i=i 



where the last equality follows from [26]. For t > 0, taking into account the definition of 
1{( T ) given at Lemma I3TT1 we obtain that Y17=i v of % ''(£) = ^> w hich contradicts the fact 
that W 1 coincides with the space spanned by 

{(/(0° /(O"'): l<j<d}. 

So we have that the Malliavin matrix Q t is invertible for any t G (0, T], as we wished to 
prove. □ 

4.3. Smoothness of the density of the solution. In order to avoid lengthy lists of 
hypothesis on our coefficients, we focus in this section on the example of the weighted 
delay treated at Section 13.31 As usual in the stochastic analysis context, we study the 
smoothness of the density of the random variable under consideration by bounding the 
L~ p moments of its Malliavin matrix. Towards this aim, it will be useful to produce an 
equation solved by the Malliavin derivative of the solution y t of equation (}4"Qj) . This is 
contained in the following Lemma: 

Lemma 4.5. Under the conditions of Proposition \3.1fy let y be the solution to equa- 
tion fjU\ ). Assume furthermore that B is a fBm with Hurst parameter H > H Q , where H Q 
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is defined at Remark \3.15[ Then y t G D 1,p for any p > 1, and $t(r) := V r y t is the unique 
solution to the following equation: 



n d 



Mr) = [U {T \y)}t + V t (r), where Vf'(r) = E E / (([W (T) (yFf ($ m ^))) s dB l s , 

m=l 1=1 ^ r 

(45) 

with the additional constraint $t(r) = for all < t < r < T. 



Proof. The equation followed by Vy is a direct consequence of relation (1431) and Proposi- 
tion [3TT21 The fact that y t G when H > H stems now from Proposition 13.141 

□ 

Now we are able to state the second main result of this section, for which we need an 
additional notation: for two a non-negative matrices M, N G M nxn , we write M > N 
when the matrix M — N is non-negative. 

Theorem 4.6. Let f, a, v and B as in Lemma \47o*[ Assume that a has bounded deriva- 
tives of any order and that 

o{rh)oW > eId R", for all 771, 772 G R n . (46) 
Then, for t G (0, T], y t has a C°° -density. 

Proof. The proof follows closely the lines of [E3, Theorem 3.5], which is classical in the 
Malliavin calculus setting, and we shall thus proceed without giving too many details. 
Nevertheless, we shall divide our proof in two steps. 

Step 1: Let Q t be the Malliavin matrix of y t , defined by (}4"4l . The standard conditions 
to verify in order to get a C°° density are: (i) y t G D°°, and (ii) [det(Qt)] -1 G LP for all 
p > 1. Condition (i) is obtained by iterating the derivatives of y, similarly to what is 
done in [26], so that we will focus on point (ii). 

In order to check that [det(Qt)]^ 1 G L p , we bound P(\ [det (Qt)] I" 1 > p) for p large 
enough, and invoke the fact that 



P (IfdetCQt)]]- 1 >p)<P \Q t % ^Idn») . 



In the sequel of the proof, we will evaluate the right hand side of the above inequality. 

Step 2: In order to bound Q t from below, the basic idea is to use decomposition (T451) 
for the Malliavin derivative of y. In this decomposition, the term [W^ T ^(y)] t is bounded 
deterministically from below under the non-degeneracy condition (I4"6j) . while V is a highly 
fluctuating quantity, since it is given by a stochastic integral with respect to B. 

One can formalize the previous heuristic considerations in the following way: 

Lt = \\uW(y)i m C = P<h (u {T \y)MoMl Hm ^ ■ 

Thanks to relation (T4~ll. one can show that 



c H T f s l ~ 2H f f\r-8) H -V\u-s) H -*l 2 r a - 1 l\ H - 1 l 2 (q:q u ,e l ) dudrds, 



i=i 



Js Js 
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where {ef, I = l,...,n} stands for the canonical basis of M n , and where we have set 
q s = cr(f® h y s+ e v(d6)) as in the proof of Proposition 13.121 Therefore, condition (T4"6l) 
yields, for a constant c which may change from line to line, 

L t > ce (J* s l - 2H J* J\r - s) H ^'\u - s f-^ r ^ u ^i/2 dudrd ^ Idfin 
> cet 2H ld R n. 

According to relation IjlSi) . it is now readily checked that 

Qt > ~~ WtWn WlRn. 
Thus, for any strictly positive number a, there exists a universal constant c such that 

It is now enough to observe that i?[||Vt||^] is a finite quantity for any p > 1, owing to 
Proposition 13. 141 to conclude the proof. 

□ 

Remark 4.7. As mentioned before, the restriction H > H for the smoothness of the 
density of the random variable y t is due to the continuous dependence of our coefficient 
/ on the past of the solution. Indeed, in case of a discrete delayed coefficient of the form 
cr(y t , yt-n, ■ ■ ■ , Vt-r q ), with q > 1 and r\ < ■ ■ ■ < r q < h, it can be seen that equation (1401) 
can be reduced to an ordinary differential equation driven by B. This allows to apply the 
criterions given in [1 3j , which are valid up to H = 1/2. 

In order to get convinced of this fact, consider the simplest discrete delay case, that is 
an equation of the form 

£0+ / a(y t ,y t - r )dB t , < t < T, (47) 

with r > 0. The initial condition of this process is given by £ G CJ on [— r, 0], and we 
also assume that a and B are real valued. Without loss of generality, one can assume 
that T = mr for m G N*. In this case, set y{k) = {y s+kr ; s G [0, r)}, and adopt the same 
notation for B. Then one can recast ([4"T|) as 

y t (k) =y r (k-l)+ [ a(y u (k),y u (k - 1)) dB u (k), t G [0, r], k < m - 1. (48) 



Setting now y = (y(l), . . . , y{m))\ B = (5(1), . . . , B{k)) t and defining a : R m -> M m ' m 
by 

^(l), • ■ • , = DiagOOK 1 )), ■ ■ ■ , v(v( m ))), 

we can express (1481) in a matrix form as 

yt = y + / *(yu(i),...,y«M)dB„„ te[o,r}. (49) 

This is now an ordinary equation driven by a m-dimensional fBm B. Whenever |cr(?7)| > 
e > and H > 1/2, one can apply the non-degeneracy criterion of [1 3J in order to 
see that y t posesses a smooth density for any t G (0, T]. The case of a vector valued 
original equation (l4"7l) can also be handled through cumbersome matrix notations. As far 
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as the case of a coefficient a(y t , yt- ri , • • • , Ut-r q ) is concerned, it can also be reduced to an 
equation of the form fl49l) by introducing all the quantities 

y t (ki, k 2 ,..., k r ) = Vt+j: r j=1 fc^-r^), 
where we have used the convention r = 0. 

Acknowledgment. Part of this work was done while Jorge A. Leon was visiting the Uni- 
versity Henri Poincare (Nancy) and Samy Tindel was visiting Cinvestav-IPN. Both are 
grateful for the hospitality of the respective institutions. 

References 

[1] E. Aids, J. A. Leon and D. Nualart (2001): Stochastic Stratonovich calculus for fractional Brownian 

motion with Hurst parameter less that 1/2. Taiwanese Journal of Mathematics 5, 609-632. 
[2] F. Baudoin and L. Coutin (2007): Operators associated with a stochastic differential equation driven 

by fractional Brownian motions. Stoch. Proc. Appl. 117, 550-574. 
[3] L. Coutin and Z. Qian (2002): Stochastic analysis, rough path analysis and fractional Brownian 

motions. Probab. Theory Relat. Fields 122, 108-140. 
[4] T. Cass, Z. Qian and J. Tudor (2009): Non-Linear Evolution Equations Driven by Rough Paths. 

Arxiv Preprint. 

[5] A. Deya and S. Tindel (2009): Rough Volterra equations 1: the algebraic integration setting. Stoch. 
and Dyn. 9 (3):437-477. 

[6] A. Deya and S. Tindel (2008): Rough Volterra equations 2: convolutional generalized integrals. Arxiv 
Preprint. 

[7] A. Deya, M. Gubinelli and S. Tindel (2009): Non-linear rough heat equations. Arxiv Preprint. 
[8] M. Ferrante and C. Rovira (2006): Stochastic delay differential equations driven by fractional Brow- 
nian motion with Hurst parameter H > 1/2. Bernoulli 12 (1), 85U-100. 
[9] M. Ferrante and C. Rovira (2009): Arxiv Preprint. 
[10] M. Gubinelli (2004): Controlling rough paths. J. Funct. Anal. 216, 86-140. 

[11] M. Gubinelli, A. Lejay, S. Tindel (2006): Young integrals and SPDEs. Potential Anal. 25(4), 307-326. 
[12] M. Gubinelli, S. Tindel (2008): Rough evolution equation. Preprint larXiv : 0803 . 05521 [math. PR], 
to appear in Ann. Prob. 

[13] Y. Hu, D. Nualart (2007): Differential equations driven by Holder continuous functions of order 

greater than 1/2. Stochastic analysis and applications, 349-413, Abel Symp. 2, Springer-Berlin. 
[14] S. Kusuoka (1982): The non-linear transformation of Gaussian measures on Banach space and 

absolutele continuity (I). J. Fac. Sci. Univ. Tokyo IA 29, 567-597. 
[15] S. Kusuoka, D. Stroock (1984): Applications of the Malliavin calculus. I. Stochastic analysis 

(Katata/Kyoto, 1982), 271-306, North-Holland Math. Library 32, North-Holland. 
[16] M. Ledoux, Z. Qian, T. Zhang (2002): Large deviations and support theorem for diffusion processes 

via rough paths. Stoch. Proc. Appl. 102, 265-283. 
[17] A. Lejay (2003): An Introduction to Rough Paths. Seminaire de probabilites 37, vol. 1832 of Lecture 

Notes in Mathematics, 1-59. 
[18] T. Lyons and Z. Qian (2002): System control and rough paths. Oxford University Press. 
[19] R. Marty (2005): Asymptotic behavior of differential equations driven by periodic and random pro- 
cesses with slowly decaying correlations. ESAIM: Probability and Statistics 9, 165-184. 
[20] S.-E. A. Mohammed (1998): Stochastic differential systems with memory: theory, examples and 

applications. In Stochastic Analysis and Related Topics VI (L. Decreusefond, J. Gjerde, B. Qksendal 

and A.S. Ustiinel, eds), Birkhauser, Boston, 1-77. 
[21] A. Neuenkirch, I. Nourdin, A. Rofiler, S. Tindel (2009): Trees and asymptotic developments for 

fractional diffusion processes. Ann. Inst. Henri Poincare Probab. Stat. 45, no. 1, 157-174. 
[22] A. Neuenkirch, I. Nourdin, S. Tindel (2008): Delay equations driven by rough paths. Electron. J. 

Probab. 13, no. 67, 2031-2068. 
[23] D. Nualart (2006): Malliavin Calculus and Related Topics. Springer- Verlag. 



30 



JORGE A. LEON AND SAMY TINDEL 



[24] D. Nualart (2003): Stochastic integration with respect to the fractional Brownian motion and appli- 
cations. Contemporary Mathematics 336, 3-39. 
[25] D. Nualart and A. Ra§canu (2002): Differential equations driven by fractional Brownian motion. 

Collect. Math. 53 (1), 55-81. 
[26] D. Nualart and B. Saussereau (2008): Malliavin calculus for stochastic differential equations driven 

by a fractional Brownian motion. Stochastic Process. Appl. To appear. 
[27] V. Pipiras and M.S. Taqqu (2000): Integration questions related to fractional Brownian motion. 

Probab. Theory Relat. Fields 118 (2), 251-291. 
[28] L. Quer-Sardanyons and S. Tindel (2007): The 1-d wave equation driven by a fractional Brownian 

sheet. Stochastic Process. Appl. 117(10), 1448-1472. 
[29] F. Russo and P. Vallois (1993): Forward, backward and symmetric stochastic integration. Probab. 

Theory Relat. Fields 97, 403-421. 
[30] J. Teichmann (2009): Another approach to some rough and stochastic partial differential esquations. 

Arxiv Preprint. 

[31] S. Tindel, I. Torrecilla (2009): Fractional differential systems for H > 1/4. Arxiv Preprint. 
[32] L.C. Young (1996): An inequality of the Holder type, connected with Stieljes integration. Acta Math. 
67, 251-282. 

[33] M. Zahle (1998): Integration with respect to fractal functions and stochastic calculus I. Probab. 
Theory Relat. Fields 111, 333-374. 

Jorge A. Leon: Depto. de Control Automatico, CINVESTAV-IPN, Apartado Postal 14-740, 07000 
Mexico, D.F., Mexico. Email: jleon@ctrl.cinvestav.mx 

Samy Tindel: Institut Elie Cartan Nancy, B.P. 239, 54506 VandoelJuvre-les-Nancy Cedex, France. 
Email: tindel@iecn.u-nancy.fr 



